Photocurrent in conjugated polymers 
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Nonlinear photocurrent carriers in conjugated polymers, such as polarons, bipolarons and solitons, 
are considered at low photon energies where a tunnelling process is necessary. We show that polarons 
usually dominate the photocurrent I due to a novel electric field assisted tunnelling for which 

In/ E-'^'^. For near degenerate polymers an electric field E which exceeds the confinement 

potential and frequencies above twice the soliton energy, soliton tunnelling is favored. Photocurrent 
data can then be used to identify the remarkable phenomenon of soliton conduction. 
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Photoexcitation of charge carriers in conjugated poly- 
mers is of considerable interest for determining types of 
charge carriers and for probing novel conduction mecha- 
nisms. Considerable theoretical attention has been given 
to polyacetylene (PA) [1-4] which has degenerate ground 
states (i.e. by inverting the dimerization sign) predict- 
ing photocurrent carried by charged solitons. Experi- 
mentally, however, significant photocurrent has been ob- 
served in non-degenerate conjugated polymers such as 
polydiacetylene (PDA) [5,6]. 

In a non-degenerate polymer the dimerization pattern 
in the ground state is unique, being fixed by the poly- 
mer structure. The charge carriers are expected to be [4] 
either polarons or doubly charged bipolarons; the con- 
tinuum model for polarons and bipolarons has been con- 
sidered in detail [7-9]. Reversing the sign of the dimer- 
ization leads to a local minimum (at least for near de- 
generate systems) which is higher by an energy of a per 
unit length. A soliton is an excitation which interpolates 
between states with opposite signs of dimerization; hence 
generating a soliton pair leads to an excess energy linear 
in the separation R of the solitons, i.e. a confinement po- 
tential aR. In contrast, polarons or bipolarons are local 
deformations of the ground state, hence a polaron pair 
has no confinement potential. We note that interchain 
coupling acts as a confinement potential, hence even PA 
with degenerate minima should be considered as a near- 
degenerate polymer. 

Some properties of PDA serve us as a prime example 
for illustrating phenomena in a non-degenerate polymer. 
PDA has chains of the form RC - C = C ~ CR =\n 
where R is one of a number of side groups with which 
the monomer can be synthesized. Band structure cal- 
culations [10] show that the states near the gap are es- 
sentially Hz orbitals. If these orbitals were removed the 
underlying structure would be [—RC — C = C — Ci?— ]„, 
i.e. 3 bonds per carbon. If a is the mean spacing of car- 
bons, the electron spectrum in an extended zone scheme 
has gaps at ±n7r/4a with n = 1,2,3 since the unit cell 
has 4 carbons. Filling in one electron per Carbon for the 
■Kz orbitals up to wavevectors ±7r/2a leads therefore to a 
gap at the Fermi level, a so called "extrinsic gap" Ag. In 



addition, the tt^ orbitals tend to dimerize, i.e. their over- 
lap between nearest carbons is alternating. This dimer- 
ization increases the gap at the Fermi level so that the 
gain in electron energies overcomes the cost in the lattice 
distortion. For PDA, the usual acetylenic ground state 
is then formed with the above triple bond. It is possible, 
however, depending on the side groups R, on temperature 
or on external strain, that the ground state will favor the 
opposite sign for this dimerization, leading to the buta- 
trienic form [—RC = C = C = Ci?— ]„. In fact, some 
data supports a PDA butatrienic form produced in a dif- 
ferent photopolymcrization procedure, which transforms 
into an acetylenic form when annealed above ~ S50K 
[11]. Solitons in PDA interpolate between these states, 
i.e. [... ~ RC ^ C ^ C ^ CR- RC -C = CCR = ...]. 
The central carbon C has only three bonds, i.e. it ac- 
quires a localized orbital with intragap energy. If a phase 
transition in PDA is indeed posible, then at the critical 
temperature a = and solitons can lead to photocurrent 
[12]. As we show below, even a small a can demonstrate 
the remarkable phenomenon of soliton conductance. 

Experimental data on PDA [5] with currents measured 
down to 10~^^A shows a steep photocurrent edge at pho- 
ton energy of 0.8eV which is well below either the band 
gap of 2.4eV, or the exciton level at 2eV. Further data 
on PDA with a different side group [6] has shown that 
thermal annealing allows for measurable photocurrents 
at 2.2eV, with an exciton level at 2.3eV; the current sen- 
sitivity in this case was about 10~^^A. 

In the present work we study photocurrent due to tun- 
nelling into a variety of charge carriers - polarons, bipo- 
larons and solitons. We show that at low photon ener- 
gies tiuj < 2Ep polarons usually dominate, where Ep is 
the polaron formation energy. The electric field allows 
for a novel tunnelling route where weak lattice deforma- 
tions are formed at large separation Rp such that upon 
charge transfer shallow polarons are formed; the energy 
gain eERp compensates for the missing formation en- 
ergy 2Ep — TiLu. This process is likely to be effective 
also in semiconductors with higher dimensionality. We 
also show that bipolarons can be directly photoproduced, 
though with a smaller probability. Finally we show that 
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the external electric field E can overcome the confine- 
ment potential when eE > a and allows for charged soli- 
ton creation. We find that for 2Es < hu) < 2Ep, and 
above the threshold field a/e soliton formation is favored, 
where 2Es is the soliton pair energy (in addition to the 
aR term). This process can serve as a clear identification 
of soliton conduction, adding to intriguing PA data [4]. 

We consider first photocurrent due to polarons at 
hw < 2Ep. The tunnelling process is parameterized by 
two local deformations of the dimcrization pattern sep- 
arated by distance R, each deformation yields a bound 
state of size I/k. These localized states are pulled from 
the valence and conduction band edges at ±Ao, respec- 
tively. We assume that direct electron-electron (e-e) in- 
teractions are embedded in the parameters such as Aq, 
Ep and k; where needed we address e-c interactions ex- 
plicitly. For the purpose of illustration, we show the elec- 
tron level structure in Fig. 1 for the electron-lattice (e-1) 
model in the absence of e-c interactions. 

The initial state potential Vo{k,R) describes the de- 
formed lattice with the electron occupation following adi- 
abatically that of the ground state. The c;xcited state 
corresponds to a dipole allowed transition of one elec- 
tron into a state with energy higher by 2uio, e.g. for 
K = we have 2luo = 2Ao. The excited state can relax 
adiabatically into a two polaron state P^P~ with energy 
minimized at some k = Kp, i.e. the polaron energy Ep{K) 
relaxes to Ep = Ep{Kp). In the e-1 model the initial state 
of each deformation has the lower level (Fig. la) doubly 
occupied while the upper level is empty; the excited state 
has one more electron in the upper state {P~) or one less 
in the lower state {P'^). 




(a) (b) 

FIG. 1. Two localized states with intragap energies lizLJo 
(continuum states starting at ±Ao, Ao > too, are not shown). 
The dots represent occupation of these levels in the e-1 model 
for (a) a separated polaron pair P~ , (b) a separated bipo- 
laron pair B'^B~ . 

The condition for having a localized state with well de- 
fined charge is (i) R^ 1/k and (u) that E does not mix 
the state with continuum states, i.e. the energy shift of 
the localized state « eE/n is small relative to the exci- 
tation energy into the continuum, i.e. eE/n <C Ag — Wq- 
We claim that for small k the latter is Ag — wq ~ Aqk^^^ 
where ^ = vp/^o and vp is the Fermi velocity in the ab- 
sence of an energy gap. For the noninteracting case this 
is seen by continuing analytically the momentum k ^ in 
in the spectrum Wq = -^/Aq — v'pK? « Ao(l — ^k^^^). 
In the final paragraph before the conclusions we show 



that to first order in Coulomb interactions the form is 
maintained. This is rather surprising since the Hartree 
term by itself gives a ~ k term [4]. The perturbation 
parameter being g^k/Ao implies that higher order terms 
lead to higher powers of n. 

The condition eE/n <C Aq — becomes then k 
where 

K„ = (eEC/Ao)'/VC- (1) 

Since typically Aq « leV, ^ « Inm and E < lO^V/ cm 
we have cE^/Aq ^ 10~^ which serves as a small parame- 
ter. This condition for localized charges is extremely im- 
portant - it allows the excited state potential to gain an 
electric field energy Vex{K, R) = 2Ep{K) — eER, hence the 
energy gain eER facilitates tunnelling even if hiv < 2Ep. 
The initial state potential has neutral components, hence 
Vo{n,R) = 2Ep{i,)-2oja(n). 

Perturbation theory in the exciting photons ~ e*"* 
[3,13] in the adiabatic limit (Tiw, 2Ao large compared 
with a typical phonon energy Tiujph) shows that the 
electron transition occurs at k,R such that Vo(k, i?) = 
Vex{i^,R) — TvjJ- The tunnelling barrier changes then 
from Vo(k, i?) to Vex{t^TR) — T'-t^' at the crossing point. 
The dynamics is dominated by the ion kinetic energy; 
the space dependent dimerization pattern A(x; fc, R) = 
Ap(K, X — R/2) + Ap(K, X + R/2), where Ap(K, x) is the 
single polaron pattern, leads to the kinetic energy 

f A^ 

/ dx— 2- = Mi{k)R^/2 + M2{k)k^/{2k^) (2) 

J nXvpuipf^ 

where dot is a time derivative and A is the electron- 
phonon coupling. Mi(k) is the mass for polaron motion 
while M2{n) is the mass for the polaron size 1/k modu- 
lations. Note that Ap{K,x) is symmetric in x hence Or 
terms are antisymmetric and the coefficient of a cross 
term Rk vanishes. 

Consider first a tunnelling trajectory such that initially 
K increases (at small i? < ^) until crossing into the ex- 
cited potential occurs. The matrix element for the transi- 
tion involves the photon field and a wavefunction overlap 
which is large for i? < ^. The tunnelling continues then 
on the excited state with increasing i?, reducing the ex- 
cited state energy with the —eER term until it vanishes 
at Ri = {2Ep-huj)/eE (Fig. 2a). The condition i?i > C 
is satisfied if 2Ep — hui is not too small, e.g. > O.leV and 
with ^ « Inm one needs E <C lO^V/ cm. Most of the tra- 
jectory is then in the regime R^ ^ justifying the use of 
the form Vex {k, R) = 2Ep — eER and the use of effective 
masses. The tunnelling rate along the R trajectory under 
the barrier 2Ep—eER—Jiu) (curve p in Fig. 2a) is given by 
a WKB form exp{-2 J^" [Mp{2Ep - eER - fiuj)]^/'^dR} 
where Mp = Mi{Kp) is the polaron mass; hence 

T, ~ e-^(2ir,-M^/= . (3) 
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This exponent involves two large parameters; first an 
electric field term, Ao/eES, > 10-^, and second an adia- 
batic parameter ■^/Mpv'^/Ao ^ Ao/tOph ~ 10. Hence Fi 
is extremely small. 
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FIG. 2. Energy barriers for (a) tunnelling along R for po- 
larons (curve p with 2Ep — Tiw — eER) and for bipolarons 
(curve 6 with 2^;, — huj — 2eER), and (b) tunnelling along 
k: the initial deformed ground state energy (curve g with 
2Ep{k) — 2ciJo{k)) crosses into the polaxon state (curve p with 
2Ep(K) — eER — TiLj) at the minimal k value of Km- 

To find the optimal tunnelling trajectory in the i?, k 
plane and avoid the strong suppression of tunnelling 
along R (Eq. 3) we consider the opposite extreme where 
K increases from at a large fixed R. This corresponds to 
nucleating two lattice deformations at a large distance R 
with initial energy 2Ep{k) — 2ujo{k) which crosses the ex- 
cited state energy at some k (Fig. 2b). The price to pay 
in this process is that the matrix elements involve the 
overlap of localized wavefunctions separated by a large 
R, i.e. InF w —2kR. As discussed above, the minimal 
K that is consistent with localized charges is hence 
we choose to nucleate the deformations at a distance Rp 
such that the crossing in Fig. (2b) is at Km, i-e. 

2Ep{Km) - 2uJo{Km) = 2Ep{Km) - hu - cERp . (4) 

Since Km is small Rp w (2 Aq — fiu) /eE, hence the tun- 
nelling rate is 

Tp ~ e-*(^-5s^^^^^ (5) 

Fp involves also the tunnelling across small k values of 
Fig. 2b which contributes a small correction within the 
exponent in Eq. (5). We note that if Coulomb energies 
would give an excitation energy ~ k (e.g. in semicon- 
ductors with higher dimensionality) then Km ^ E^^^ and 
In/ ^ (Ao/ei?^)^/^. We also note that choosing R < Rp 
for Fig. 2 implies a higher k so that the above KmRp is 
the minimal value. 

The result for Fp is remarkable - it reduces the large 
Ao / eE£^ factor in Eq. (3) to its 2/3 power and also avoids 
the large adiabatic parameter. Furthermore, this process 
is much more efficient than the Franz-Keldysh process 
[14] which is the photon assisted k = process, i.e. tun- 
nelling between extended states of the valence and con- 
duction bands with InF (1 - ?iaj/2Aof/'^Ao/eE^. 



We conclude that the process leading to Eq. (5) is the 
most efficient one for generating polarons. This is a most 
remarkable process - the system prepares localized states 
at a large separation (typically Rp > lOOnm) so as to ac- 
commodate an eventual transfer of charge. 

We consider next tunnelling into bipolarons. The elec- 
tric field breaks inversion symmetry and allows charge 
transfer between polarons, leading to doubly charged 
bipolarons. The excitaion potential is Vb{K, R) = 
2E\j{k) — 2eER where 2Eb is the creation energy of two 
bipolarons; in the e-1 model 2Ei,{k) = 2Ep{K) + 2loo (Fig. 
lb). If tunnelling along R had dominated, then direct 
tunnelling into bipolarons would have dominated at low 
frequencies. To see this, note that Vi){k, R) — hu> crosses 
at ii2 = (2i^6 — ^) /2 {Eb here is the minimized Eb{K)) 
which is smaller than Ri (Fig. 2a) if fioj < AEj, — 2Ep, 
or h(jj < 0.4Ao in the e-1 model. Tunnelling directly into 
bipolarons eliminates the more costly tunnelling range 
between i?2 and i?i at the expense of an e~2Ki?2 factor 
for the second charge transfer. Hence, bipolarons would 
dominate at low Ttw. 

Within the much more efficient process of tunnelling 
along K at a large R bipolarons are more efficiently 
generated when R = Rb is chosen such that the bipo- 
laron curve 2Eb{K) — 2eER — h,uj crosses the initial en- 
ergy 2Ep{k) — 2w(y{K) at the minimal Km-, hence Rb ~ 
(2Ao — \huj)/eE. The polaron excitation energy is then 
below the bipolaron one, but since —eER is inefficient 
at K < Km the polaron curve would also cross near Km- 
The transfer of two charges involves now a e"^**^ factor, 
hence F^, <C Fj, with 

r,^e-«(i-^»^5^'\ (6) 

We consider next photocurrent of charged solitons 
which has been studied in detail for the degenerate PA 
case [f 4]. In the non-degenerate case, or for PA with in- 
terchain coupling, the electric field must exceed a thresh- 
old value to overcome the confinement potential. Soli- 
tons being topological objects must be generated by R 
tunnelling - there is no k parameter or a local defor- 
mation which can generate a single soliton from the 
ground state. The excitation energy of a soliton pair 
S+S- is then Vs{R) = 2Es + aR - eER allowing tun- 
nelling only for eE > a. The tunnelling terminates at 
Rg = {2Es — hu))/{eE— a) which is assumed to be large, 
Rs ^ ^. The tunnelling rate into solitons with mass Mg, 
for hu) < 2Es, eE > a, is then 



The formation of S^S~ leaves a metastable state which 
can then spontaneously form a pair of solitons with- 
out a photon. Fg is extremely small (comparable to 
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Eq. 3) unless Tilo is close to 2Es- Soliton pliotocurrent 
may then dominate the polaron one only if the range 
Eg < Tujj < Ep exists. In the e-1 model with weak a one 
has Eg ~ O.6A0 < Ep ^ O.9A0 (Es is defined so that a 
soliton pair energy separated by i? » ^ is 2Eg + aR) . 

We consider therefore the range 2Es Kfiu) < 2Ep. The 
polaron photocurrcnt is still given by Eq. (5) while soli- 
ton tunnelling can occur at i? < ^ where a and E are 
ineffective. The result is then just as in the a = sys- 
tem [3] InF « —Aq/ujq, with weak E dependence. We 
conclude that soliton photocurrent dominates over that 
of polarons when 2Es < hu) < 2Ep and E > a/e. This 
is achievable if a/e < lO^V/ cm or a = AAeAo/nXvp [9] 
yields 7 = Ag/AAo ^ 0.1. Hence varying E changes Eq. 
{5) at E > a/e into a larger and weakly E dependent 
tunneling rate, exhibiting a sensitive tool for identifying 
the intriguing soliton conduction. 

Finally, we present a somewhat technical section, eval- 
uating the effect of Coulomb interactions to first order 
on the excitation energy of a polaron, e.g. P+, lead- 
ing to the important criteria for Km in Eq. (1). The 
e-1 system with (non-uniform) electron transfer between 
nearest neighbors has particle-hole symmetry, i.e. each 
eigenstate (f>{n) at site n with energy E allows an eigen- 
state (— with energy —E. Completeness rela- 
tion for the polaron P+ state (right side of Fig. la) 
yields for the charge density of each spin Pi{n) = ^, 
Pli^) — \ ~ 'Pbi''^) where 4)b{n) is the lower energy lo- 
calized eigenstate (Fig. 1). A Hubbard interaction, to 
first order is X^„(Pt(") ~ i)(Pl("') ~ ^) = 0; similarly all 
interactions with range of even sites vanish. 

Consider next nearest neighbor interaction TLy with 
coupling V. Its average yields direct and exchange terms, 

(Hy) = iy^{(p(n)-i)(p(n + l)-i) 

n 

-|5 + 0^(n).^6(n+l)|2-|5|^} (8) 

where p{n) = p|(n) + Px("-) and S = + 1) 

sums on all occupied continuum states (per spin) with 
energy up to — Aq. The 2nd term above is the exchange 
for 1 spin while the 3rd is for J, spin. The summation can 
be done in a continuum limit [9] where x = na and , 
are wavefunctions on the even and odd sites, respectively. 
For small k {1/k is the range of the localized state) we 
obtain from [9] 

s = - (1 - H>^)h{x)m{x) + oiK-") (9) 

and the contribution to the polaron energy becomes 

{-Hv), = Vr + Va ^: 7 ^° (1 - -7^^) + 0{^) (10) 

where Vi/N = -Vo?[{Ae - Aq)/2-k\vfY is the correc- 
tion to the ground state energy per site. Note that the 
direct term 5 V /fc ('^•) ~ 0{k) cancels with with 



the exchange term [Ref. [4] presents just the direct term 
which, as shown here, is insufHcient] . 

The polaron P+ excited state involves an electron 
transfer from the top of the valence band with an ex- 
tended wavefunction 0„(n) into the lower localized state 
(Fig. la) which becomes doubly occupied. The direct 
term becomes | X^n '/'^ (") 'Z'? ('^+ 1 ) ~* for an extended 
system, the t exchange is the same as in Eq. (8) while the 
J, exchange becomes \S+(f)l{n)(t)b{n+l)-(j)%{n)(t>v{n+l)\'^. 
Using [9] / dxf*{x)gy{x) = ^ + 0{k'^) we obtain that the 
excited polaron energy is the same as that of the polaron 
to order k^. We have also checked this conclusion within 
the continuum model with c-e interactions. Hence the 
excitation energy is 0{k^), the same as in the e-1 model. 

In conclusion, we have shown that polarons usually 
dominate the photocurrent when hoj < 2Ep with the tun- 
nelling rate from Eq. (5), i.e. In / ~ —E~^^^. However, 
for 2Es < fkv < 2Ep and E > a/e solitons dominate 
the photocurrent. This is achievable for near-degenerate 
polymers, i.e. a small a. Furthermore, the intriguing 
possibility that PDA has a phase transition at which the 
confinement parameter a vanishes [11] can be sensitively 
tested by the field dependence of photocurrent data. 
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